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ABSTRACT 

We present the results of a spatio-kinematic study of the planetary nebula Menzel 1 using 
spectro-photometric mapping and a 3-D photoionization code. We create several 2-D emission line 
images from our long-slit spectra, and use these to derive the line fluxes for 15 lines, the Ha/H/3 
extinction map, and the [SII] line ratio density map of the nebula. We use our photoionization 
code constrained by these data to derive the three-dimensional nebular structure and ionizing star 
parameters of Menzel 1 by simultaneously fitting the integrated line intensities, the density map, 
and the observed morphologies in several lines, as well as the velocity structure. Using theoretical 
evolutionary tracks of intermediate and low mass stars, we derive a mass for the central star of 
0.63±0.05M Q . We also derive a distance of 1050±150pc to Menzel 1. 

Subject headings: Planetary Nebula - Interstellar medium 



1. Introduction 

The Planetary Nebula (PN) Menzel 1 (Mz 1) or 
G322.4-02.6 (15 h 34 m 16 s . 7 -59° 08' 59" 2000.0) is 
a bright object with a bipolar morphology and a 
prominent central ring of enhanced emission. Ha 
and [OIII] narrow band images of Mz 1 have been 
published by Schwarz, Corradi, & Melnick (1992). 
Being bright has not resulted in Mz 1 being well- 
studied; only a few papers have been dedicated to 
the object. 

H2 emission has been detected from Mz 1 by 
Webster, Payne, Storey, Dopita (1988) and the 
morphology in this molecular line is similar to that 
of our optical image, shown in Fig. 1. 

A detailed kinematical study of Mz 1 was done 
by Marston, et al. (1998) using high resolution 
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echelle spectroscopy. Using the intensity ratio of 
the [SII] lines they computed densities of 1700 
cm~ 3 for the main ring and 400 cm -3 for the bipo- 
lar part. They show that the gas in the main part 
of the nebula is expanding at 23 km s , and that 
the velocity structure is consistent with a cylindri- 
cal model with expansion velocities proportional 
to the radial distance from the center. They also 
show that the dynamical age of the ring in Mzl 
is of the order of 7000 yrs and estimate a mass of 
about 0.5M Q for the nebula using dynamical ar- 
guments and their computed distance of 2.0 kpc. 

Distances have been determined by several au- 
thors such as Van de Steene, & Zijlstra (1995) with 
2.53 kpc, Cahn, Kaler, & Stanghelini (1992) with 
2.28 kpc and Zhang (1995) with 2.85 kpc, all using 
statistical methods and based on the same data. 
Using the same formalism as Van de Steene, & Zi- 
jlstra (1995) and a new value for the radius of the 
nebula, Marston et al. (1998) obtained a distance 
of 2.0 kpc with a claimed uncertainty in this es- 
timate of about 30%. Acker et al. (1992) lists 10 
distances of which 8 are statistical and 2 individ- 



ual determinations. Given that all the above au- 
thors use a filling factor, e = 0.5, any such method 
has an inherent and large uncertainty in the dis- 
tances they determine. We compute 2.0±0.5kpc 
(adjusted standard deviation) from all 14 litera- 
ture distances without any weighting factor. Note 
that each individual distance can be in error by a 
large factor (see the appendix), and the standard 
deviation computed from the literature values is 
not very robust. 

An effective temperature of 139 kK and a lumi- 
nosity of 147 L Q were determined by Stanghcllini, 
Corradi, Schwarz (1993) with a distance of 1.8 kpc. 
Mz 1 is clearly under-luminous for a typical PN. 

In this work we present observations of, and 
a 3-D photoionization model for Mz 1 , and de- 
rive the 3-D structure of the PN constrained by 
observed fluxes and morphologies in many emis- 
sion lines, using the same method as Monteiro, 
Schwarz, Gruenwald, & Heathcote (2004) applied 
to NGC 6369. By determining the 3-D structure of 
nebulae, the large uncertainty involved in all clas- 
sical statistical distance determination methods is 
eliminated. Assuming an arbitrary filling factor, 
constant ionized mass or diameter, mass-radius re- 
lationship etc. is not needed here: we determine 
what the structure and ionized mass are, and can 
therefore derive distances to much greater accu- 
racy than has been previously possible. 

In summary, we obtain the 3-D spatial struc- 
ture of the nebula along with its chemical com- 
position, ionizing source temperature, luminosity, 
and mass, as well as an independent distance, in a 
self-consistent manner. In §2 we discuss the obser- 
vations and briefly explain the basic reduction pro- 
cedures, including our image reconstruction tech- 
nique used to obtain the emission line intensity 
images. In §3 the results obtained from these im- 
ages are discussed: the pixel by pixel reddening 
correction of the images, the integrated line fluxes, 
and the computed temperature and density maps. 
In §4 we present the model results generated by 
our 3-D photoionization code, and we discuss the 
derived quantities. In §5 we give our conclusions, 
and explain in detail how our method works in the 
appendix. 



2. Observations 

2.1. Observations and data reduction 

We show our image of Mz 1 in Fig. 1 taken in 
the light of the [SII]671.7nm line through a fil- 
ter with A c = 671. 8nm and FWHM = 2.6 nm. The 
300 s exposure was taken with the CCD camera 
attached to the CTIO 0.9m telescope on the 9th 
of April 2002. The plate-scale is 0.4"/pix on the 
2kx2k TEK chip, and the seeing during the ex- 
posure was 1.1" according to the nearby seeing 
monitor. 
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Fig. 1. — Our narrowband [SII] image of Mz 1 with 
the scale and orientation indicated. We have indi- 
cated the centers of the slit positions we used by 
the short strokes on both sides of the image. 
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The main features of this [SII] image are sim- 
ilar to those found in the Ha image of Schwarz, 
Corradi, & Melnick (1992). The morphology of 
Mz 1 is complex: two faint outer lobes toward the 
NW and SE meet in the brighter central annular 
region, and there are enhanced extended bands of 
emission on the eastern and western side of the 
central annulus. Some fine structure is also seen 
throughout the surface of the nebular emission. 

The spectra were taken with the CTIO 1.5 m 
Ritchey-Chretien telescope with the RC Spectro- 
graph on the nights of 13 and 14th June, 2002. 
We used a grating with 600 1/mm blazed at 600 nm 
giving a spectral resolution of 0.65nm/pixel and 
a plate scale of 1.3"/pixel with a slit width of 
4". The spectral coverage obtained with this con- 
figuration was approximately 450 nm to 700 nm. 
For details of the instrument and telescope see 
http://www.ctio.noao.edu and click on "Optical 
Spectrographs" , then on "1.5m RC spectrograph" . 

By taking exposures at several parallel long-slit 
positions across the nebula, we obtained line inten- 
sity profiles for each slit. These profiles were then 
combined to create emission line images of the 
nebula with a spatial resolution of about 4" x 4", in 
a way similar to radio mapping. Due to a techni- 
cal problem the data from one of the observed slit 
positions had to be discarded. We computed an 
average of the two adjacent slit spectra as repre- 
senting this position (shown in Fig. 1). The added 
uncertainties introduced by this procedure are dis- 
cussed below. 

The individual slits were observed and re- 
duced using standard procedures for long-slit spec- 
troscopy, using IRAF reduction packages. 

A fine correction for slit misalignment was made 
using the Ha and H/3 profiles for each exposure. 
Using IDL, the images were rc-dimensioned to 100 
times their original size. The normalized Ha and 
H/3 profiles were then matched and the final result 
re-dimensioned to original values. This procedure 
yields the precise alignment necessary for calcu- 
lation of diagnostic line intensity ratios. Minor 
shifts of the order of one pixel can introduce con- 
siderable errors in line ratios, if this method is not 
applied. 

It is also important to note that the 11 slit 
positions observed do not fully cover the nebula, 
leaving out the faintest outer parts of the bipolar 



lobes; see Fig. 1. Below we estimate this lost flux 
and show that it is a small fraction of the total 
flux. Note, however, that in our calculations, we 
take this "lost flux" into account by matching our 
model output fluxes to the observed area, but the 
model does provide the complete nebular struc- 
ture. 

3. Observational results 

Images were created for all 15 lines detected 
with a signal-to-noise ratio above 5. In Fig. 2 the 
constructed images for the most important emis- 
sion lines are shown. These images have been cor- 
rected for reddening as described in the following 
section. The corresponding signal-to-noise images 
were used to obtain the total line intensities with 
fractional errors, given in Table 1. 
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3.1. Reddening correction and total fluxes 



From our long-slit spectra we reconstruct emis- 
sion line images using the method described 
by Monteiro, Schwarz, Gruenwald, & Heathcote 
(2004). These images for each line were corrected 
for reddening using the Ha/H/3 ratio map shown 
in Fig. 3. The logarithmic correction constant 
was calculated pixel by pixel using the theoretical 
value of Ha/H/3=2.87 from Osterbrock (1989) and 
the reddening curve of Seaton (1979). 

We investigated the effect of differential atmo- 
spheric refraction on this ratio map. From the 
airmasses of our observed positions and the values 
given by Filippcnko (1982) we computed a correc- 
tion which we applied to our data. All slit posi- 
tions were observed at airmasses below 1.5 except 
one outer position in a faint part of the nebula, 
which had an airmass of 1.8. Since we used a 
4"wide slit and the object is extended, the effect 
was small, but not negligible in the steep gradients 
near the bright ring structure. The average error 
due to this effect is about 3% in the high signal to 
noise (S/N) areas and about 25% in the low S/N 
areas. The nett effect on the final calculated rela- 
tive total fluxes is about 0.3% for strong lines and 
7% for weak ones, well within the other observed 
uncertainties. 
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Fig. 2.— Reddening corrected emission line images for Mz 1. Only the 10 images with the highest S/N are 
shown. 
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Fig. 3.— The Mz 1 Ha/H/3 ratio map with con- 
tour overlay of Ha, where the darker shading rep- 
resents stronger emission. Only the brighter parts 
of the nebula are shown as the image is cut for 
S/N values lower than 10. 



7 



The final calculated line fluxes relative to H/3 
and their corresponding la errors are shown in 
Table 1 . All total fluxes were obtained by integrat- 
ing the reddening corrected images pixel by pixel. 
The value we obtained for the total reddening cor- 
rected H/3 flux is F Hfj = 2.6 xlO-^erg cm- 2 s- 1 . 
We estimate from the relative area of the neb- 
ula that falls outside the slit images, and the 
fact that for several lines the flux comes mainly 
from the central, bright part of Mzl, that the 
lost flux in H/3 is about 3% of the total, and 
the average lost flux in all lines about 5%. Note 
that Acker et al. (1992) lists an uncorrected H/3 
flux of 4.9x 10~ 12 erg cmT 2 s _1 , taken from Perek 
(1971), which is close to our uncorrected value of 
5.8 x 10 _12 er<? cmT 2 s _1 . Our flux is obtained dig- 
itally by integrating our spectrophotometry pixel 
by pixel over the whole nebula, while Perek used 
photoelectric aperture photometry during non- 
photometric nights, with all its associated calibra- 
tion problems, and warns us in his article that cir- 
rus affected many of the measurements, explaining 
why his flux is lower than ours which was taken 
during a photometric night. Perek also listed a 
previous measure of the H/3 flux of -11.26 which is 
equal to our uncorrected flux. 

Using these images we estimated the effect of 
adopting the average for the discarded slit posi- 
tion. Since the nebula shows considerable sym- 
metry, we compared the northern half of the re- 
constructed images with the southern half. By 
comparing the total fluxes from these halves we 
determined that the error on the calculated to- 
tal fluxes were about 10% for low ionization lines 
and about 2% for the other ones. This difference 
is due to the fact that the discarded position is 
close to the bright ring of the nebula, where the 
low ionization lines are relatively strong. The ef- 
fect is relatively small and this added uncertainty 
is taken into account in the values presented in 
Table 1. 

3.2. Gas density and temperature 



Heathcote (2004). The density map is shown in 
Fig. 4. 

For the temperature maps, the correction for 
slit misalignment was carried out as done for the 
Ha/H/3 maps discussed in section 2.1. We com- 
pute the temperature map for Mz 1 in two ways: 
using our density map we obtain the temperature 
map shown in Fig. 5 (A) and in (B) we show the 
difference between this map and the one calculated 
assuming constant density. The maximum differ- 
ence between these two maps is less than 100 K. 
The images are clipped for data values with S/N 
lower than 10 for visual clarity. 



We calculated density and temperature maps 
from the reddening corrected maps of the [SII] and 
[Nil] lines respectively. The expressions relating 
the line intensity ratios to the gas density and tem- 
perature are the ones published by McCall (1984) 
and those used in the IRAF temden package. De- 
tails are as in Montciro, Schwarz, Gruenwald, & 
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Table 1: Line fluxes relative to H/3 



Line A(nm) 


Flux 


Dered. Flux 


Error (%) 


HcIIA468.6 


0.11 


0.11 


7 


[OIILU495.9 


2.35 


2.28 


6 


[OIII]A500.7 


7.11 


6.83 


6 


[NII]A575.5 


0.10 


0.07 


23 


HeIA587.6 


0.25 


0.16 


9 


[OI]A630.0 


0.33 


0.21 


8 


[SIII]A631.1 


0.03 


0.02 


20 


[OI]A636.3 


0.12 


0.08 


14 


[NII]A654.9 


2.96 


1.81 


12 


Ha 


4.68 


2.88 


6 


[NII]A658.4 


9.18 


5.57 


12 


HeIA667.8 


0.08 


0.05 


15 


[SII]A671.7 


0.82 


0.49 


12 


[SII]A673.1 


0.76 


0.45 
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Fig. 4. — Density map obtained from the 
[SII]671. 7,673. lnm line ratio. 
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Fig. 5. — (A) Temperature map obtained from 
the de-rcddcncd [N II] line ratio. The scale on 
the right is in degrees centigrade; (B) difference 
between this map and the one obtained assuming 
constant density. Note that the maximum differ- 
ence is less than 100 K. North is up, and east is to 
the left. 
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4. Photoionization Models 



The photoionization code applied here has 
been described in detail by Gruenwald, Viegas, 
& Broguiere (1997). It uses a cube divided into 
cells, each having a given density, allowing arbi- 
trary density distributions to be studied. Typical 
models runs use cubes of 70 to 100 cells on a side. 
The input parameters are the ionization source 
parameters (luminosity, spectrum, and tempera- 
ture), elemental abundances, density distribution, 
and the distance to the object. The conditions are 
assumed uniform within each cell for which the 
code calculates the temperature and ionic frac- 
tions. These values are used to obtain emission 
line emissivities for each cell. 

The final data cube can be spatially oriented 
and projected in order to reproduce the observed 
morphology. The orientation on the sky of the 3- 
D nebular structure is thus determined. The line 
intensities and other relevant quantities are then 
obtained after projection onto the line of sight. 

The structure we obtained for Mz 1 is shown in 
Fig. 6 as it is oriented relative to the observer. We 
also show a cut along the major axis of symmetry 
in Fig. 7, indicating the density values. The 3-D 
structure of Mz 1 was constrained by our density 
map to be an open structure, not a closed shell of 
any type. We therefore adopt an hour-glass shape 
with a density gradient from the equator to pole. 
We added random density fluctuations to better 
fit the line fluxes, especially H/3. The rotation an- 
gles relative to the x, y, and z axes respectively, 
are 0°,10°, and 40°, with the symmetry axis of 
the main hour-glass structure being x. The model 
resolution of 100 3 cells was limited by our 4 GB of 
computer memory, and code execution time. 
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Fig. 6. — Isocontour plot of the 3-D density struc- 
ture determined for Mz 1 using model calculations. 
The inclination angle of the polar axis of the struc- 
ture is 40°. 
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The ionizing spectrum we used for the central 
star (CS) was a blackbody modified by the He and 
H atomic absorption edges at 54.4eV and 13.6eV. 
The addition of these edges was necessary to be 
able fit the [OIII] and Hell line intensities simulta- 
neously, the Hell line being particularly sensitive 
to the depth of the absorption edge, and therefore 
providing a good central star temperature con- 
straint. These absorption edges are also present 
in the more sophisticated theoretical atmospheric 
models presented by Rauch, Dectjcn, Drcizler & 
Werner (2000), which are similar to the ones we 
used here for similar stars. The only feature from 
Rauch's spectra that we do not incorporate is the 
added absorption due to other lines, but this has 
no significant effect on our model. The effective 
temperature and luminosity of our adopted crude 
stellar spectrum are given in Table 2. 

4.1. Model Results 

We present here the main results obtained from 
the photoionization model constrained by the ob- 
servational data. The total line intensities are 
given in Table 2, as well as the fitted abundances 
and ionizing star parameters. Projected line im- 
ages for four important transitions are shown in 
Fig. 8. 

The model image size is fitted to the ob- 
served one for the line [NII]658.4nm, as well as 
the absolute HP flux, giving a final distance of 
1050±150pc. Fig. 9 shows the observed [Nil] im- 
age with the corresponding model image contours 
overlaid for the obtained distance. 
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Fig. 8. — Projected line images obtained with the photoionization model. 
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Fig. 9. — Contour image of model [Nil] image 
overplotted on observed image for our best dis- 
tance of d=1050±150pc. 



Using the model cube of temperatures and ion- 
ization structure we also calculated the position- 
velocity (PV) diagram for four different slit posi- 
tions using a spherically symmetrical velocity field 
given by V — a(\r\ jrmax) x r/|rj . r max is half 
the side of our model data cube (here we used 
Ymax = 8 10 17 cm), and a is the maximum veloc- 
ity reached within the data cube (here we used 
45km/s). We computed the PV diagram for the 
[NII]658.4nm line in the same positions as those 
observed by Marston et al.(1998). The projected 
velocities in the x direction are obtained using: 



with: 

AV(v, x, y, z) = V x (x, y, z)-v 



(1) 



£(x,V,z) = \/v t 2 h (x,y,z) + V 1 



(2) 



(3) 



V th (x,y,z) = y / 2kT e (x,y,z)/Am H 



(4) 



where e\(x,y, z) is the emissivity of a cell at 
{x,y, z), Vt is the turbulent velocity (taken to be 
2km/s) and Vth(x,y, z) the thermal velocity of an 
atom of atomic mass A, the local electron temper- 
ature is T e (x,y, z), and v is the plotting interval 
for the velocity, where in this case v goes from -80 
to +80km/s. irh is the mass of an H atom. The 
PV diagram obtained is shown in Fig. 10. 



Table 2: Observed and model line fluxes and model 
central star parameters. 





Observed 


Model 


T* (K) 




120kK 


L-j, 1 Lq 




164 


Density 


100-1400 


100-1400 


He/H 




1.14 x lO" 1 


C/H 




3.3 x 10~ 4 


N/H 




2.2 x 10~ 4 


O/H 




4.7 x 10~ 4 


Ne/H 




3.5 x 10~ 4 


S/H 




1.1 x 10~ 5 


\og(H(3) 


-10.6 


-10.6 


[NeIII]386.8 a 


5.64 


4.8 


HcII468.6 


0.014 


0.013 


[OIII]500.7 


6.83 


6.68 


HeI587.6 


0.16 


0.17 


[OI]630.0 


0.21 


0.26 


[NII]658.4 


5.6 


5.3 


[SII]671.7 


0.49 


0.49 


[SII]673.1 


0.45 


0.46 



a) Value obtained by Kaler, Shaw, Browning 
(1997) 
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Fig. 10. — Model PV diagram obtained for the 
four slit positions as observed by Marston, et al. 
(1998) in the light of the [NII]658.4nm line. 



19 



Comparing our model derived PV diagrams 
with the observed ones in Figure 3 of Marston 
et al.(1998) we see a remarkable agreement. The 
main features are reproduced, confirming our con- 
fidence in the 3-D structure we adopted for Mz 1. 

5. Discussion and conclusions 

We present spectrophotometric maps of Mzl, 
giving spatially resolved information for many 
emission lines and precise total fluxes for the ob- 
served part (>95% of the total flux) of the neb- 
ula. The images produced with this technique 
were used to determine properties of the nebula 
using the Ha/H/3 ratio map for the de-reddening, 
the [SII] line ratios for the density, and the [Nil] 
line map for the temperature. 

The Ha/H/3 map shows little structure. The 
prominent ring feature of the nebula does not show 
significant differences in extinction when com- 
pared to the other regions. This may imply that 
most of the reddening is due to the foreground, 
and not intrinsic to the nebula. 

We also show the density map obtained from 
the observations indicating the presence of a dense 
waist ring and a bipolar structure of lower den- 
sity. Based on this map we propose for Mz 1 a 
3-D hour-glass structure with a waist whose den- 
sity decreases smoothly from the equator to the 
poles. 

Using a photoionization code and the proposed 
structure we obtained a complete 3-D physical 
model for Mzl. The fitted model line intensities 
show excellent agreement -well within the com- 
puted errors- with the observed values. The ob- 
tained distance of d=1050±150pc falls near the 
lower end of the error range determined from 14 
literature values for the distance. This makes 
sense as the angular radius of 25" used in the sta- 
tistical methods is smaller than the true size of the 
nebula of 38" as measured by Schwarz, Corradi, 
& Melnick (1992), and our distance is also smaller. 

The luminosity of 164 ± 25 L and temperature 
of 120 ± 16 kK we have determined for the central 
star are within their errors equal to those found by 
Stanghellini, Corradi, Schwarz (1993), taking into 
account the difference in distance and total flux. 
Our errors have been conservatively estimated to 
be similar to the cumulative observational errors 
we computed and used to constrain the model. 



This luminosity is low for a PN CS, and indi- 
cates that the star is evolved and well down its 
Schonberner track, as confirmed by its position in 
Fig. 11. 
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Fig. 11. — Comparison of our model temperature 
and luminosity obtained for the central star with 
model tracks calculated by Blocker (1995). 
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We determine the mass of the CS to be 
0.63±0.05M© by fitting to the theoretical evolu- 
tionary tracks of Blocker (1995) , using the errors 
estimated from our observational data and propa- 
gated to the CS luminosity and temperature. Due 
to the fact that in this part of the HR diagram, 
the tracks lie close together, we cannot say with 
any precision what the precursor mass of this star 
was, except to place a lower limit on this mass 
of about 1 M Q using the added track from Vassil- 
iadis, & Wood (1994) in Fig. 11. The likely value 
for the CS precursor mass is about 3M , based 
on the most probable CS mass of 0.63 M©. From 
our model output we compute the ionized nebular 
mass to be 0.14 ±0.03 M . 

An independent check is the time scale of ex- 
pansion of Mzl. From Marston, et al. (1998) and 
scaling to our distance we get an expansion time of 
3500 yrs for the ring which can be longer or shorter 
depending on what the history of the expansion 
has been -decreasing velocity due to energy con- 
servation in a wind blown cavity would shorten 
the time and increasing velocity due to expan- 
sion into a medium with decreasing density would 
lengthen it-so a range of about 2000-5000 yrs. For 
the outer parts of Mzl this time would be about 
4700-12000yrs as the material has had to travel 
2.4 times further away from the CS (38" instead 
of 16"). From our luminosity range and comput- 
ing Blocker track ages we obtain 4500 to 10000 yrs, 
quite compatible with the outer nebular expansion 
time, thus confirming our distance determination. 

Since the sum of the ionized nebular mass and 
the CS mass is 0.77 M©, we expect that there pos- 
sibly is more than 2 M© of neutral matter near 
Mzl. Clearly, for the lower CS precursor mass 
limit of 1 M© , the neutral mass estimate is a fac- 
tor of ten smaller. In any case, there is likely more 
neutral than ionized mass in the system. Note, 
however, that Huggins, Bachiller, Cox, Forveille 
(1996) estimate a molecular mass of 0.0086 M© 
and an ionized mass of 0.027 M©. They used a 
radius for Mzl of 12.6" instead of the 38" optical 
radius, so naively assuming that the mass scales as 
the third power of the radius would increase their 
ionized mass by a factor of up to 9 to 0.25 M©, 
which is much closer to our computed mass of 
0.14M©. 

We have put together the Spectral Energy Dis- 
tribution (SED) from 0.36 to 100 fiia for Mz 1, us- 



ing data from the literature listed in SIMBAD. 
The plot of XF(X) is shown in Fig. 12. Note that 
the 100 /im IRAS point is an upper limit. The 
double peaked distribution is typical for a PN. 
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Fig. 12.— The SED for Mz 1 between 0.44 /an and 
100 yum using flux values taken from the literature. 
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Integrating our SED yields a luminosity of 
about 60 L Q which includes the usual correction 
proposed by Myers et al. (1987) of a factor of 1.5. 
This implies that about 60% of the UV radiation 
from the CS escapes from the nebula. Given that 
the structure we found is open toward the poles, 
and is likely clumpy, this is a reasonable value. 
The fact that we estimate that there is between 
0.2 and 2 M Q of neutral matter around Mz 1, some 
of this may have, in fact, become ionized due to 
this escaping radiation. Searching for a large, low 
surface brightness halo around Mz 1 may therefore 
be useful. 

Perhaps most important is the fact that our 
simple method can be applied to any spatially re- 
solved emission line nebula, yielding the complete 
3-D structure and an accurate distance. The time 
consuming multiple position long slit observations 
can be dispensed with since Integral Field Spec- 
trographs are now becoming available at major 
observatories, making the method not just sim- 
ple but also efficient. See the detailed explanation 
of our method in the appendix. 

6. Appendix: Determining distances us- 
ing 3-D photoionization models 

Classical distance determination methods are 
statistical or individual in nature and all assume 
constancy of one or more physical parameters of 
the PNe. Gurzadyan (1997) provides an excellent 
review of distance determination methods and we 
refer to his book for details on other methods than 
the "Astrophysical Method" (AM) which is the 
one we use here. 

The AM uses the fact that the electron density, 
H/3 flux, and angular extent are observable quan- 
tities of PNe -treated as spherically symmetrical 
objects- and that they are related to the distance 
of the nebula by: 

d = 2.4.10 25 F(H/3)/[n 2 .6» 3 .e] 

where F(H/3) is the H/3 flux in ergs/cm 2 /s, n 2 
is the electron density in cm~ 3 , 9 is the observed 
angular extent of the nebula, and e is the so-called 
filling factor which is the fraction of the nebular 
volume which is emitting i.e. it contains ionized 
gas. 

It is therefore in principle possible, if all the 
above parameters are known, to determine an ac- 



curate distance to any resolved, spherical nebula 
or HII region. In practice, however, the method 
has severe limitations due to the necessity of re- 
strictive assumptions such as spherical symme- 
try Taking each parameter in turn, we show that 
a typical distance determination has an error of 
about a factor of 3 or more for any given nebula. 
There are nebulae for which individual distances 
with a range of a factor of 100 have been published 
(M2-9 extremes are: 50 to 5200 pc; Schwarz, As- 
pin, Corradi, & Reipurth (1997))! 

H/3 fluxes are typically measured in a small 
aperture centered on the PN, and an attempt is 
made to estimate the flux from the entire neb- 
ula by extrapolation. Typical errors are large due 
to incomplete knowledge of the size and bright- 
ness structure of the nebulae. The presence of the 
central star in the aperture increases the uncer- 
tainty in this measurement. Estimated errors are 
in the range 1.5 for smooth, regularly shaped neb- 
ulae without much structure to 3 or more for more 
typical PNe. 

Electron densities, n, are usually determined 
from the ratio of the pair of forbidden sulfur lines 
at 671.7 and 673. lnm. Since spectra are also typ- 
ically taken with small apertures centered on the 
PN or in the best cases with a long slit but for 
only one position across the nebula, again the un- 
certainty in the average value of n is large, say, a 
factor of 1.5-2. Note that n appears in the equa- 
tion to the second power, increasing the effect of 
this error. 

The angular extent, 9, of a nebula can be mea- 
sured but depends on the passband or emission 
line used to make the observations. A PN can 
have a diameter in the light of the [Nil] line that is 
twice that measured in H/3 or [OIII] , again produc- 
ing uncertainties in the distance determination. 
Smaller nebulae have a correspondingly larger un- 
certainty associated with this measurement. 9 ap- 
pears in the equation to the third power, increas- 
(fr)g the effect of the uncertainty in this parameter. 
An error of about a factor of 1.5 - 3 is again typical. 
In any case, the use of the radius assumes spherical 
symmetry, which for most PNe is far from realis- 
tic. 

The most difficult to determine and therefore 
least known parameter is the filling factor, e. Usu- 
ally, it is taken to be 0.5 in statistical methods 
using large samples of nebulae, but the true value 
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can be anywhere between and 1, with typical es- 
timates varying between 0.2 and 0.9, i.e. a factor 
of about 5 uncertainty. 

Propagating all these factors, it is clear that 
determining a distance to a PN is a very uncertain 
business. 

Our method differs fundamentally from the 
classical methods, in which one or more of the 
above parameters were assumed to be constant 
or have some simple relationship. The use of 
a 3-D structure to model the nebula eliminates 
the need for assumptions on the quantities n,6, 
among others, as the structure can be modified to 
include small scale variations in density such as 
clumps, filaments and others as well as large scale 
variations (hour-glass shapes for example). More 
importantly, this procedure removes entirely the 
need to specify a filling factor e, as the large and 
small scale density variations are all well defined 
in the 3-D structure. We also determine all of the 
observables with high precision from either long 
slit or Integral Field Unit (IFU) spectra across the 
nebula. Our photoionization model is then con- 
strained by the quantities derived from these de- 
tailed spectra, and spectral images. We constrain 
simultaneously with: several line images, several 
line fluxes, complete projected density map, and 
the velocity structure, obtaining the best overall 
fit by adjusting the central star luminosity, spec- 
tral distribution, temperature, average chemical 
abundances, and the distance, also obtaining the 
complete and detailed 3-D structure of the nebula. 
All the above mentioned parameters are therefore 
known much more precisely, and the distance de- 
termination is correspondingly better. Typically 
we can compute the distance to about 10-20% de- 
pending mainly on the observational errors. 

6.1. Model fitting procedure 

The details of the numerical photoionization 
code are given in the appendix of Gruenwald, Vie- 
gas, & Broguiere (1997). In summary, the gaseous 
region, with the radiation source in the center, is 
divided into cubic cells, for each of which the phys- 
ical conditions are, by definition, homogeneous. In 
each cell thermal and ionization equilibrium is as- 
sumed in order to obtain the physical conditions. 
The radiative transfer problem is calculated with 
the "on the spot" approximation to save comput- 
ing time. 



The input parameters are the elemental abun- 
dances, the gas density distribution, the shape and 
intensity of the central ionizing radiation spectrum 
(temperature, H and He absorption edges, and lu- 
minosity in the case of a star) and the distance to 
the object. The code then provides the physical 
conditions in each point of the nebula, i.e., ionic 
fractional abundances, electronic temperature and 
density, as well as the emission line luminosities of 
each cell. 

This output is then used to produce projected 
images and total fluxes for a given number of emis- 
sion lines and a set of (x,y,z) orientation angles. 
This output can then be tailored to match given 
observational configurations such as single long 
slits, or multiple slits (as is the case for Mzl). 
From the projected images we can also construct 
projected diagnostic maps, such as density and 
temperature maps. 

The model generated or simulated "observa- 
tions" discussed above are then compared to the 
actual observational data obtained (in this case, 
the spectro-photometric mapping of Mzl). We 
compare total line intensities and check for dis- 
crepancies. If one or more of the intensities are out 
of the range of the observational errors, we proceed 
to fine tune input parameters that have influence 
on the given line. For example, we take the Hj3 
total fluxes from the model and observations and 
compare them. The HP line is mainly dependent 
on the star luminosity and the 3-D structure of the 
gas, so we adjust these input parameters accord- 
ingly. In this case it is important to realize that 
the 3-D structure is actually defined by the den- 
sity in each cell and a physical size for the object 
which is dependent on the distance. So in fact we 
are dealing with two input parameters when we 
consider the H/3 flux. The same type of compari- 
son is made for other lines such as [OIII]500.7nm 
which is an important coolant, HeII468.6nm which 
depends mainly on the CS temperature, among 
others. We also compare the model projected im- 
ages and diagnostic maps to those obtained from 
the observations. 

After fitting all the model constraints to their 
respective observational counterparts and adjust- 
ing the input parameters of the code accordingly, 
we calculate a new model. The same procedure 
discussed above is then repeated until we reach a 
satisfactory agreement between model and obser- 
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vations for all line images, fluxes, diagnostic maps, 
etc. 

Notice that, after this iterative procedure, we 
obtain model fitted values for the input param- 
eters that are self-consistently determined; they 
are the ionizing star characteristics, gas chemical 
abundances, density, structure, and distance. 

One of the main differences of this procedure 
when compared to previous model calculations in 
the literature, is that we use the distance as a fit- 
ting parameter for the model. This is possible be- 
cause we now have a way of producing projected 
images from the 3D model results and can there- 
fore compare them directly with observed ones as 
well as the observed total fluxes. In other words, 
we do not use a fixed distance for our model calcu- 
lations and vary only star and gas parameters (lu- 
minosity and temperature of the star, abundances 
and densities of the gas) . The other important ad- 
vantage of using the 3-D structure for the gas is 
the possibility of eliminating the need for a "filling 
factor" . This has major implications on the model 
parameters that can be determined, especially on 
the distance. 

Since we use a 3-D structure that is consistent 
with observed images, position-velocity diagrams, 
diagnostic ratios (such as density maps) we are not 
making any assumptions about filling factors, ion- 
ized masses or physical sizes. All these parameters 
are determined self-consistently in the model. 

The uniqueness of the solution obtained by this 
procedure can be argued of course. It is imme- 
diately clear that within the observational uncer- 
tainties there are an infinite number of solutions 
that can fit the observations. In other words, the 
observations determine the quality of the final pa- 
rameters. In fact, if we estimate the goodness of 
fit of our model fit by quadratically summing all 
uncertainties and dividing by the number of ob- 
servables minus the number of degrees of freedom 
(input parameters to the model), we get an er- 
ror of about 20% in the case of Mzl. This is the 
uncertainty adopted for our results. The precise 
determination of fitting errors is extremely com- 
plex and given the nature of the 3-D code neither 
practical nor useful. 

For the above reasons our distance determina- 
tions arc fundamentally different from, and much 
more precise than classically found distances. 
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